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Abstract: In this paper, Landau theory for phase transitions is shown to 
be a useful approach for quantal system such as atomic nucleus. A detailed 
analysis of critical exponents of ground state quantum phase transition 
between 1/(5) and 0(6) limits of interacting boson model is presented. Our 
results suggest a similarity between these two frameworks and a second 
order phase transition between these limits based on a discontinuity in the 
heat capacity. 
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1. INTRODUCTION 

Studying the behavior of nuclear matter under extreme conditions of 
temperature and density, including possible phase transitions, is one of the 
most interesting subjects in recent years. Drastic changes in the properties 
of physical systems are called phase transitions which these properties have 
been characterized by order parameters. Phase transitions occur as some 
of parameters, i.e. control parameters, which have constrained system, 
are varied. Temperature-governed phase transitions in which the control 
parameter is temperature, T, have been known for many years [1], Landau 
theory of phase transitions [2-3] was formulated in the late 1930’s as an 
attempt to develop a general method of analysis for various types of phase 
transitions in condensed matter physics and especially in crystals. It relies on 
two basic conditions, namely on (a) the assumption that the free energy is an 
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analytic function of an order parameter and on (b) the fact that the expression 
for the free energy must obey the symmetries of the system. Condition (a) is 
further strengthened by expressing the free energy as a Taylor series in the 
order parameter. 

For fluid systems, as we become close to the critical point, some of the quantities 
of system are related to the temperature as f(T) « (T 7’) ,; for some exponents 
of (3 . Thesimilar behaviors may be seen not as a function of temperature but as 
a function of some other quantities of system, e.g. f(x)«(xf. These exponents, 

(3 , are called critical exponents and naturally defined as lim —l—l ; [4], Some 

ln(±x) 

basic critical exponents in thermodynamics have been employed to describe the 
evolution of considered systems near the critical points [5-6]. 

2. QUANTUM PHASE TRANSITION IN THE INTERACTING 
BOSON MODEL (IBM) 

In nuclear physics, quantum phase transitions, sometimes called zero 
temperature or ground-state phase transitions can be studied most easily 
with using algebraic techniques that associate with a specific mathematical 
symmetry with different nuclear shapes. Interacting Boson Model (IBM) as 
the most popular algebraic model in description of nuclear structures was 
proposed in 1975 by lac hello and Arima to describe the collective excitations 
of atomic nuclei [7-10]. In this model, nucleons in an even-even isotope are 
divided into an inert core and an even number of valence particles. These 
particles are then considered as coupled into two kinds of bosons that may 
carry either a total angular momentum 0 or 2, and are respectively called the 
s- and 7-bosons. The bilinear operator that may be formed with s- and J-boson 
creation and annihilation operators close into the (7(6) algebra whose three 
possible subgroup chain match with the (7(5), SU(3) and SO(6) solution of 
the Bohr Hamiltonian, i.e. respectively with spherical, axially deformed and 
'-'(-unstable shapes. It is of great interest to be able to describe the evolution of 
considered systems near the critical points. Let’s consider a general form of 
IBM Hamiltonian as [7]. 


H = e d h d + k 0 Q x=0 .Q x=0 +k l p\p + k 2 L.L + k 3 T 3 J 3 + kj 4 .T A 


( 1 ) 


where h d {=d\ d) is the d boson number operator and Q x °, i.e. 

(s f x d + dd x i) (2) explores the quadrupole interaction. Also, other terms of 
Hamiltonian are 
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This general Hamiltonian can describe three dynamical symmetry limits 
with different values of constants, i.e. e 0 , X and k. s. We must consider a 
transitional Hamiltonian to describe the critical exponents at the critical 
point of phase transition. To this aim, we propose the following schematic 
Hamiltonians for U( 5) —0(6) transition [11,21] 


H t ra„ S =^h d +^P\P, (3) 

Where h d (= d' . d) and P f = -T. T) are introduced in 

Eqs(l-2) and we have considered e d = x and k 1 = (l-x)/(N-l). The 
U( 5) limit of IBM is recovered via x = l and x = 0 reproduces the 0(6) 
limit. It means one can describe a continuous, e.g. second-order shape 
phase transition by changing x between these two limits. On the other 
hand, classical limit of transitional Hamiltonian, Eq.(3), is obtained by 
considering its expectation value in the coherent state [12-14] 

\N,a m ) = (s + '£a m d m ) N |0), (4) 


Where |o) is the boson vacuum state, ,v' and of are the creation operators 
of s and d bosons, respectively and a m can be related to deformation collective 
parameters, a 0 = (3 cos 7, a ±1 =0 and a ±2 = /? sin 7/V2 . The energy surface 
which follows from expectation value of transitional Hamiltonian in the 
coherent state, Eq.(4), is given by 


E(N,/3) = N[x 


( 3 2 

1 + /3 2 


1 — x 

4 


j 1 P 2 ?] 
I + P 2 


(5) 


Critical point of considered transitional region have obtained via 
d 2 E(N,/3)/d/3 2 \^ 0 = 0 [15] condition which gives x critical = 0.5 in this case. We 
show the dependence of energy surface on the order parameter, f3 , above 
and below of the critical point of phase transition, x , in Figure 1. In 
phase transition from, U( 5) i.e. spherical limit, to 0(6), namely, 7-unstable 
limit, one sees that, the evolution of energy surface goes from a pureto [i 1 
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Figure 1: Energy surfaceswhich are determined via transitional Hamiltonian, Eq.(l). 
Different panels describe the dependence of energy surfaces on the order parameter, (3, 
above and below of phase transition point, x 

A r critical 


a combination of f3 2 and /? 4 that has a deformed minimum. At the critical 
point of this transition, energy surface is a pure ,3 4 . These results interpret 
that d 2 E(N,/3)/d(3 2 1 o = 0 condition corresponds approximately to a “very 


flat energy surface”, similar to what have happened for the E(5) critical 
point [16], i.e. critical point of 17(5) *->0(6) transitional region. The quartic 
dependence of the energy at the critical point and the predictions for spectra 
and B( E2) transition rates for the energy surfaces /3 4 ,/3 s ,/3 8 which are called 
respectively, the E(5)-(3* , E(5)- /3 6 , E(5)-(3*, have been given in Refs. 
[22-23]. Some comparison of the work performed in the present letter with 
the variational method used in literature ([49], for example) is in place. 
The main difference between the two methods is that in Ref. [22,24] a trial 
potential (the Davidson potential) containing a parameter( (3 0 ) is chosen and 
subsequently the rate of change of the physical quantity (the rate of change 

E(L) 

of the energy ratios, R = ——' ) is maximized with respect to this parameter, 

E( 2) 

thus determining the parameter value and the value of the physical quantity 
(the energy ratios) while in the present letter a method based on Catastrophe 
Theory in combination with a coherent state analysis has been employed to 
generate the energy surfaces and the evolution of the potential energy surfaces 
with respect to (3 are considered. Also, The analysis carried out in Refs. 
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[25-26] concentrated on the spectroscopic properties of the Ru isotopes and 
on the application of nuclear supersymmetry for a simultaneous description 
of these nuclei. Here we have considered the geometrical behavior of 
U( 5) «-> SO (6) transitional Hamiltonian. 

The typical behavior of the order parameter, p , at a phase transition is 
shown in Figure 2. Here B is small and close to x . and we assume that 
energy surfaces can be expandedaround /? = 0, 


Critical Exponents 
of Quantum 
Phase Transition 
Between U(5) and 
0(6) Limits of 
Interacting Boson 
Model 


1 _ y 

E(N,p)=[N( - ) + N(2x-l)p 2 + — (—3 lx+l9)P 4 +....]. (6) 

4 12 


Or can be rewritten in the form 

E(N,a,b,c m ,B) = a + bp 2 +c/3 4 , (7) 


The behavior of, near the critical point is determined by the signs of the 
coefficients b,c. The coefficients b,c which are functions of x, are written as 
functions of the dimensionless quantity x = {x — x critical ) / x critical , 


2 ml order 




**,cnicm 05 


Figure 2: Typical behavior of order parameter, (3, at a second order phase transition. 


a = N 


l — x 


b = Nx, C = — (7-3 lie) 
24 


( 8 ) 


where x critical = 0.5, b is represented as a series in odd powers of Jc so that 
its sign would change at c s = 0 . On the other hand, the stable systems must 
have same signs on both sides of x critical = 1/2 which this requirement is equal 

to c > 0 . This means is represented only as 2_n ■ 

24 
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The condition for stability is that, E(N,a,b,c-,(3) must be a minimum as a 
function of f3 . From Eq. (7), this condition may be expressed as where terms 
above ,3 4 are presumed negligible near x critical = 0.5 [17]. For x critical < 0.5 
, only real root is /? = 0 ; on the other hand, for x critical > 0.5 , the root j3 = 0 
correspond to a local maximum, a^d therefore not to equilibrium. The other 

two roots are found to be j3 2 =- x . Consequently, our analysis predicts, 

the equilibrium order parameter nenr the critical point should depend on the 
x as 


/3=Constant x (— x) 2 , (10) 

which means, critical exponent for order parameter is 1/2. The behavior 
of (3 0 (x) is depicted in Figure 3, which is in perfect agreement with general 
predictions derived in Ref. [2] . 

On the other hand, a very sensitive probe of phase transitional 
behavior is the second derivatives of the ground state energy (per boson) with 
respect to the control parameters [18]. 



Figure 3: The dependence of equilibrium deformation, (3 (x),(a) and the specific 
heat, C(x), (b) of the ground stateabove and below of phase transition point, x . ( . for 
second order phase transition. 


C( A,. 




e 0 ( A )» 


( 11 ) 


(all X.s with j i are kept constant). In the above discussed thermodynamic 
analogy e 0 (A) is replaced by the equilibrium value of the thermodynamic 
potential F 0 (P,T) ■ In our descriptions, by use of Eq. (7), ground-state 
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energies are a,a -for /3 = 0 and (3 ^ 0 respectively. From Eq. (11) the 

4c 

specific heats are 


C + (f3g =0) = 0 for x< 0.5 ; C + (/3 0 * 0) = y Nforx> 0.5 (12) 


These results propose any dependence of C on x either above or below of 
x „i,icai = 0-5 and therefore, the values for the specific heat exponents are both 
zero. Also, this result suggests a discontinuity in the heat capacity in the 
phase transition point which in the agreement by Landau’s theory. We have 
represented the behavior of specific heat in Figure 3, which one can find that 
it has a jump at the critical point. 

The classification of phase transitions that we follow in this paper and that 
is followed traditionally in the IBM is the Ehrenfest classification [17,19]. In 
Ehrenfest classification, first order phase transitions appear when there exist a 
discontinuity in the first derivate of the energy with respect to the control parameter. 
Second order phase transitions appear when the second derivative of the energy 
with respect to the control parameter displays a discontinuity. It can be seen from 
Figure 4, that first derive of the energy surface has a king a.tx crilical . This corresponds 
to a second order phase transition, as the second derivate is discontinuous. 

In order to identify other critical exponents, according to the Landau 
theory, by use of Eq.(7), the potential energy surface becomes as [4,20]. 
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E <P> iE ip 



Figure 4: Variation of energy surface,E((3) , and its first derivative respect to order 
parameter, dE I d/3 , which are presented above and below of phase transition point, 

^ critical' 


E((3) = a-h/3 + bl3 2 +c/3* + ■■■ (13) 

Where, h/3 , represents the contribution of intensive parameter, h, for 
points off the h = 0 coexistence curve. The equilibrium equation of state is 
^/d$ xh = 0 which cause to (for any small h) 
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h = 2 b/3 + 4c 3 


(14) 


On the other hand, it reduces to its former representation for h = 0 . The 
susceptibility may be found as it introduced in Ref. [4,20] , namely, 


dh 

d/3 


2Nx + — N/3 2 
2 


(15) 


For x < 5 which we have (3 = 0 and consequently we get X*' = 2/V.i, which 
gives the critical exponent equal to 1. For x < 0.5 with h = 0 , Eq.(13) gives 
12 

(3 2 =—-x and therefore x(' = 4N(—x) or the critical exponent equal to 1. 
Along the critical isotherm, i.e. in the phase transition point, namely x = 0 


7 

and /t =-A/3 which this means, critical exponent is equal to 3. 
6 


summarize the values of the critical exponents. 

Our results, i.e. behavior of order parameter about critical point, 
discontinuity of the second order derivative of energy respect to order 
parameter, suggest a second order shape phase transition between U(5) and 
0(6) limits of IBM. Also, critical exponent and their capability to describe 
the order of quantum phase transition may be interpreted a new technique to 
explore shape phase transitions in complex systems. 


Table 1. Critical exponents of ground state quantum phase transition between U( 5) 
and 0(6) limits. 


Critical exponent 


definition 


values of the critical exponents 


Order parameter 


(3 ~ (—xY(3 =1/2 


Specific heat for * > ~ jT“C = 0 for X < X M C ~ \x\~° C = 0 

Susceptibility for x „ (7= , for X < X^X? - 1 = 1 

Critical isotherm ^ ^ h l,s 8 = 3 


3. SUMMARY AND CONCLUSION 

In this contribution, we show that, U(5) «-> 0(6) shape phase transition are closely 
related to Landau theory of phase transition and explore some of the analogies 
with thermodynamics. Also, a detailed analysis of the critical exponents of ground 
state quantum phase transition is presented. We find that, critical exponents in 
two frameworks are similar. Based on a discontinuity in the heat capacity in the 
phase transition point, we can conclude the order of the phase transition. 
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